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An analysis is made for the variabl 

laminar free convection on an isothermal vertical f^fit plate. For a num- 
ber of specific cases, solutions of the boundarylayer equations appro- 
priate to the variable property situation were carried out for gases and 
for liquid mercury. Utilizing these findings, a simple and accurate 
shorthand procedure is presented for calculating free convection heat 
transfer under variable property conditions. This calculation method is 
well established in the heat transfer field. It involves the use of re- 
sults which have been derivedr-«£pr constant property fluids, and of a set 
of rules (called reference temperatures) for extending these constant 
property results to variable property situations. For gases, the con- 
stant property heat transfer results are generalized to the variable 
property situation by replacing, P (expansion coefficient) by l/T m 
and evaluating the other properties at T r = T w - 0.38(T W - T^) . For 
liquid mercury, the generalization may be accomplished by evaluating all 
the properties (including p) at this same T p . It is worthwile noting 
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COMMITTEE 


that for these fluids, the film temperature (with 3 = l/T — for gases 
appears to serve as an adequate reference temperature for most applica 
tions. Results are also presented for boundary layer thickness and 
velocity parameters. 

NOMENCLATURE 


an 

A 1 ,A 2 

b 

c 

C P 

Cf 


F 

8 

Gr x 


Gr 


L 


h 

h 

k 

L 


Nu x 

Nu l 

P 

Pr 

q 


coefficients in the polynomial representations of fluid 
properties of liquid mercury 

dimensional constants in Sutherland ; s formulas of table 1, 
A 1 = 362° F abs, A 2 = 198.7° F abs 

plate width 

dimensional constant defined by eq. (6a) 
specific heat at constant pressure 
coefficient of friction defined by eq. (29) 
dimensionless dependent variable defined by eq. (6b) 
acceleration due to gravity 

Grashof number based on x, dimensionless, see eqs. (13), 
(13a), (13b) 

Grashof number based on L 

local heat transfer coefficient, q/(T w - T„) 
average heat transfer coefficient, Q/Lb(T w - T_) 
thermal conductivity 
plate height 

local Nusselt number, hx/k, dimensionless 
average Nusselt number, hL/k, dimensionless 
static pressure 

Prandtl number, Cpp/k, dimensionless 

local heat transfer rate per unit area of plate 
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* 




Q 

T 

t 

Tf 

u 

^ax 

v 

W 

x 

y 

y u 

3 




T\ 

e 

p 

V 

P 

T i 

* 


r L 

over- all heat transfer rate, b / q dx 

. 0 


absolute temperature 
Fahrenheit temperature 

film temperature, (T w + Tj/2; or t f = (t w + t^)/2 
velocity component in x direction 
maximum value of u across the boundary layer 
velocity component in y direction 

rate of fluid flow generated by free convection, b / pu dy 


0 

coordinate measuring distance along plate from leading edge 
coordinate measuring distance normal to plate 


distance from plate at which u mny occurs 
coefficient of thermal expansions!, - 


boundary layer thickness defined as distance from plate at which 
(T - Tj/(T W - Tj = i 


dimensionless similarity variable defined by eq. (6a) 
dimensionless temperature variable, (T - T — )/(T - T — ) 

absolute viscosity 
kinematic viscosity, p/p 
density 

shear stress at plate surface 
stream function 


Subscripts 

w denotes wall conditions 

„ denotes ambient conditions 

r denotes conditions at reference temperature 


L 
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INTRODUCTION 

The presence of a bu6yancy force is a requirement for the existence 
of a free convection flow. Ordinarily, the buoyancy arises from density 
differences which are a consequence of temperature gradients within the 
fluid- Any analytical treatment of free convection must include density 
variations at least to the extent that a buoyancy force enters the 
problem. 

A characteristic common to previous analytical studies of free con- 
vection has been the neglect of all fluid property variations, except for 
the essential density differences noted above*. Such a simplified treat- 
ment does not appear unreasonable when the temperature differences in- 
volved are small. This intuitive feeling has been corroborated in a 
formal ma nn er by Ostrach (2). For situations where there are large tem- 
perature differences, the adequacy of the results derived from the con- 
stant property analysis has been in doubt. 

An analytical treatment of the variable property problem, including 
numerical solutions, first appears to have been given in the thesis (l) 
from which this paper is taken. Some time after the appearance of Ref. 

1, a much less extensive study was described in a Russian article by 
Tanaev (3). There is no significant overlap between Ref . 3 and the pre- 
sentation here. 

The analysis is made for an isothermal vertical plate, and the flow 
is taken to be laminar. For a large number of specific cases, numerical 

It 

The fluid properties entering the problem are the thermal conduc- 
tivity, viscosity, specific heat, and density. 


solutions of the boundary layer equations appropriate to the variable 
property situation were carried out for gases. It was initially planned 
to also study the variable property problem for several liquid metals. 
However, the nature of the property variations and the exceedingly time 
consuming numerical computations forced the current study to be limited 
to one liquid metal, namely, mercury. 

Prime attention is focused on the heat transfer. Utilizing the heat 
transfer results corresponding to special cases for which numerical solu- 
tions were obtained, a simple and accurate shorthand procedure is pre- 
sented for computing heat transfer under variable property conditions in 
gases and liquid mercury. This method is well known in the heat transfer 
field. It involves the use of results derived for constant property 
fluids, and of a set of rules for extending the constant property results 
to variable property situations. These rules are commonly termed refer- 
ence temperatures. Not only will reference temperatures be derived for 
the heat transfer, but also for boundary layer thickness and velocity 
parameters. Those who sure interested primarily in results are invited j 
to pass over the section on analysis. 

ANALYSIS 

Physical model and coordinates . - The physical model and the co- 
ordinate system are portrayed in an elevation view in Fig. 1. Two phys- 
ical situations are shown which come within the scope of the analysis. 

The left hand sketch depicts the case where the wall temperature, T w , 
exceeds the ambient temperature T . Under these circumstances the free- 
convection motion is upward as shown. The right hand sketch shows the 
situation where T w is lower than the ambient temperature T — . In this 
case, the fluid flow is downward along the plate. 
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If the coordinate systems are taken as indicated, the mathematical 
distinction between the two situations vanishes when the conservation 
equations, as written later, are made dimensionless. So, separate analy- 
ses need not be made. Since it seems easier to visualize occurrences 
associated with the hot wall cas,e, i.e., T w > the analysis will be 
directed toward that situation. However, the results will be presented 
in a manner applicable to both T w > T„ and T w < T w . 

Conservation laws . - The equations expressing conservation of mass, 
momentum, and energy for steady flow in a boundary layer on a vertical 
plate are 



( 1 ) 

( 2 ) 

(3) 


/ St . dr\ d L dl\ 

p °p( u s + V ^J- 

Viscous dissipation and work against the gravity field have been 
neglected. 

The boundary conditions appropriate to the problem are 


V = 

o - 



-“1 




u = 

o , 

u = 

0 

k y s 0 






T = 

T. 

T s 

T w^ 



- 


y *» - 


(4) 


(5) 


where T w and T — are prescribed constants. 

From Eq. (3), it follows that the pressure p is a function of x 
alone (i.e., a function only of height along the plate). So, dp/dx is 
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constant across the boundary layer and may be evaluated as 

bp/bx = - p.g 

With this substitution, Eq. (2) becomes 

p ( u S + v s g(p - " p) + If) (2a) 

There appears in this equation a density difference, p m - p, which pro- 
vides the buoyancy force for the free convection motion. 

The customary approach followed in the constant property analysis 
is first to replace the density difference by a temperature difference, 
and then to assume that p, p, Cp, and k are constant. 

The present analysis continues without placing any restrictions on 
the nature of the property variations. The solution of Eq. (l) may, as 
usual, be written in terms of a stream function i|r defined by the 
relations 


p „ _ d* 




(la) 


where p w , the fluid density at the wall, is regarded as a constant. 

Then, the velocity components u and v which appear in Eqs. (2a) and 
(4) are replaced in favor of the stream function if. The result of the 
substitution is a rather complicated looking pair of simultaneous partial 
differential equations for + and T as functions of x and y. 

Rather than deal with these two formidable partial differential equations 
directly, experience suggests a method of transforming them to ordinary 
differential equations, which are easier to solve. In the usual termin- 
ology of boundary layer theory, such a transformation is called a simi- 


larity transformation. 
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Reduction to ordinary differential equations . - A new independent 
variable tj, called a similarity variable, is defined by 

: l/4 f* P 


TJ = CX 


Jo p w 


where 


c = 


e(p. ■ p w )/pw 


dy 


1/4 


4y ‘ 


w 


f , ) 


(6r) 


New dependent variables F and 0 are given by 

^ T - T. 

- JiTl ’ 4V^’ 0 ( T l) - T w •- T. 


(6b) 


The function 6 is a dimensionless temperature and F is related to the 
velocities in the following way 

- 2 tc 1/2 F', v« — : [tjF’ - 3 f] (7) 


u = 4vc x 

W 


The primes represent differentiation with respect to tj. 

U^der the transformation Eqs. (6a) and (b), the partial differential, 
equations for y and T become 


. r ~\ ? (p./p) - i 

LBL. F" + 3FF" - 2(F’) + 7 7 — , = 0 

J W? - 1 


| P|_ e>| + 3Pr. 

dTJ 

w 


' + 3Pr t/^-V 0 ' * 0 

VPw/ 


( 8 ) 


The Prandtl number is represented by Pr and the subscript w denotes 
conditions at the wall (y = 0) . The boundary conditions, Eq. (5), trans- 
form to 

F • 0 


F' = 0 
0*1 


Vij » 0 


F 1 = 0 
0=0 


> tj = - 


O) 
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It may be noted that the transformation relations, as well as the 
resulting ordinary differential equations, can be reduced directly to the 
well-known equations for the constant property fluid. Further, all solu- 
tions of Eqs. (8) subject to the boundary conditions (9) have the charac- 
teristic that u = 0 and T = along the line x = 0. 

Simplifications associated with p = pRT . - For the special case of 
fluids which obey the perfect gas law, p = pRT, there are certain simpli- 
fications which can be introduced into the analysis. Consider first the 
buoyancy force, g(p^ - p), which appears in Eq. (2a). For a perfect gas, 
this density difference can be transformed into a temperature difference 
without making any approximations. It is only necessary to notice that 
both p^ and p are to be evaluated at the same value of x, and con- 
sequently, at the same static pressure p. With this in mind, it follows 
directly that 

P„ - p = jr(T - Tj (10) 

In the constant property analysis, it has been customary to write 

P„ - P * ft?(T - Tj 

where p is the coefficient of expansion. So, it is seen that the 
buoyancy force used in the constant property analysis is precisely cor- 
rect for a perfect gas, provided that p is replace by l/T^. 

Further, in the differential Eqs. (8), the dimensionless buoyancy 


force (the last term) can be simplified to 

( P- /p) - 1 T - T„ 
(P„/P w ) - 1 “ T w - T-“ 


(10a) 
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SOj Eqs. (8) become 




(8a) 


It is worthwhile noting that for a very special sort of variable 
property gas, Eqs . (8a) reduce to the equations for a constant property 
fluid. Consider a gas having the property variations: p = pRT, pp = 

constant, pk = constant, c p = constant. Inspection of Eqs (8a) show the 
tremendous simplification afforded by the assumption that pp = constant 
and pk = constant. It follows that 

pp/PwPw = 1, Pk/Pw^w = 1 


and also that k/p = constant. Further, it is seen that the Prandtl 
number (= Cpp/k) is also a constant. Then Eqs. (8a) may be rewritten as 

F"' + 3FF" - 2(F') 2 + 0 = 0 

(8b) 

0" + 3Pr F0 ' = o 


But, these equations are precisely those for the constant property prob- 
lem. Also identical are the boundary conditions. So, from the mathe- 
matical point of view, the constant property problem is identical to that 
for the special variable property fluid: p = pRT, pk = constant, 

pp = constant, c p = constant. All solutions which have been obtained for 
the constant property differential equations become .available for this 
special variable property fluid. 

It is interesting to note that a similar finding applies in forced 
convection. In fact, a common procedure for accounting for variable :: 

■f/ • v r ■ ■ ’ r ; ’ * ■' ■ J . 1 '* 5 . '' ’ '■'* ;■ ■*. . j. v ,, i > USt 


property effects for forced convection over a flat plate is to postulate 
that the real gas is sufficiently well approximated by the special gas 
just considered. 

HEAT TRANSFER PARAMETERS AND GRASHOF NUMBERS 
Local heat transfer . - The local heat flux from the surface to the 
fluid may be calculated using Fourier ' s Law 



Introducing the dimensionless variables from Eqs. (6a) and (b), the ex- 
pression for q becomes 

1 - - *»<*, - :.)« x/ f (11) 

The derivative [de/dtQ^Q, normally abbreviated 0'(O), is found from the 
solutions of Eqs. (8). A dimensionless representation of the results is 
achieved by use of the local heat transfer coefficient and local Nusselt 
number, which are written in the usual way as 


q hx 

h = T w " T m ’ Nu x = k ( 12 ) 

Further, a generalized Grashof number which is applicable to both con- 
stant and variable property fluids is defined by 


Gr x = 


_ g* 5 (1 p» ~ Pwi/Pw) 


( 13 ) 


The absolute magnitude sign removes the necessity for separate consider- 


ation of T w > T and T w < T — . 


Using these definitions of h, Nu x , and Gr x the local heat flux 


given by Eq. (ll) becomes 
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where the second subscript on the Nusselt and Grashof numbers indicates 
the location at which k and v are evaluated. 

Overall heat transfer . - The overall heat transfer Q is found by 
integrating Eq. (ll). So, 

Q = b q dx (15) 


where b is the plate width. Introducing the following dimensionless 
groups 

— Q — hL ^ {J P - ” p w|/ p wj ,-,-x 

h = VT l m . ~ rp I) NU t S t, , G r T. = •? (16) 


leads to the following dimensionless result for the overall heat transfer 


Nu t 


, ± 0 ( 0)3 


L,w 


L ' u ” 3 V5 

The Grashof number. - The generalized Grashof number defined by 


(17) 


Eq. (13) arises naturally from the analysis of the general variable 
property fluid. The form taken by this definition for certain special 
cases is of interest. 

For instance, for the constant property analysis, Eq. (13) simpli- 
fies to 


Gr x = 


gP T w - T. 


(13a) 


This is identical to usual Grashof number definition. 
For a perfect gas, the Grashof number becomes 

g|T w - Tjx 3 


Gr x = 


(13b) 


T v‘ 


- 13 - 


DESCRIPTION OF THE GASES STUDIED 
In previous variable properties analyses for gases, it has been 
common to use idealized forms of the property variations. There have 
been a few instances where real gas properties were used. In the present 
study, both idealized and real gases have been included, the real gas 
being a close approximation to air. Table X describes the five gases to 
be considered here in the order (reading to left to right) in which they 
will be discussed. All are seen to obey the perfect gas law, p = pRT, 
Absolute temperatures are used exclusively throughout the analysis for 
gases . 


TABLE I 

Description of Gases A, B, C, D, and E 


Gas A 

Gas B 

Gas C 


GaS D 


Gas E 

p = pRT 

p = pRT 

P : “ . 

pRT 


= 

pRT 

P 

= pRT 

k-T 3 / 4 

k ~ T 2 / 3 

pk = 

const . 

k 

w 

3/2 
T ' 

T + 

k 

m 3 / 2 

T 

* T + A ± 

p ~ T 3 / 4 

[1 m 

PM * 

const . 


«v 

T3/2 
T + A x 


T 3 / 2 
~ T + Aj, 









» b 6 + b,T 

c_ = const. 

C-n = const. 

C-n = 

const . 

C T> 

SS 

const . 



P 

p 

P 


P 



P 

“ b H + 

Pr = const. 

Pr = const. 

Pr = 

const. 

Pr 

= 

const . 

Pr 

= variable 


Gases A, B, and C represent simple idealizations of real gas be- 
havior. The power law variations for k and p. are commonly used 
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approximations. The assumption of constant specific heat and Prandtl 
number is included because, for real gases, the variation of these prop- 
erties is small compared to those of k, p, and p« The Sutherland type 
formula used to describe the conductivity and viscosity variations of gas 
D is closer to reality than are the simple power laws. Since c^ and 
Pr are maintained as constants for gas D, the same Sutherland formula 
is used for k and p. Gas E is a close approximation to air. Hence, 
variations of Cp and Pr are included, and different Sutherland for- 
mulas are used for k and p (A^ft Ag). Two linear equations are used 
to represent the specific heat variation over the temperature range 
studied. 

HEAT TRANSFER RESULTS FOR GASES 

Heat transfer results for gases A through E will be presented for a 
large number of special cases. Utilizing these results, a rapid and 
accurate shorthand method will be presented for calculating heat transfer 
to gases under variable property conditions. As described in the 
INTRODUCTION, the calculation procedure involves the use of results which 
have been derived for constant property fluids, and of a rule for ex- 
tending these constant property results to variable property situations. 
This rule is commonly called a reference temperature. 

The following approach will be used here in deriving and testing a 
reference temperature rule for heat transfer. First, using the numerous 
numerical calculations for gas A, a reference temperature rule will be 
derived. Then, tests of the wider applicability of this reference tem- 
perature result will be made using the less numerous solutions for gases 
B, C , D , and E . 
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Results for gas A . - The properties of gas A are given in Table I. 

In order to compute the heat transfer, it is first necessary to 
solve the differential Eqs. (8a) subject to the boundary conditions (9). 
The gas properties appearing in Eqs. (8a) may be evaluated for gas A with 
the aid of Table I. So, it is seen that Cp/cp = 1, Pr * constant and 


P4 _ pk _ /_T_ 


tY 1 / 4 


p |i p k 
w w w w 


(18a) 


wy 


Introducing the dimensionless temperature 6 = (T - T^)/(T W - T^) , it 
follows that 

r 1 / 4 


PH 
P M- 


pk ' 
w w 


( T.\ Tj 

[ 1 - T 1 + t 

\ v/ w 


(18b) 


In light of this, it may be observed that before commencing with a solu- 
tion of Eqs. (8a) to find F and 0, it is necessary first to specify 
the value of T w /T w (as well as of the Prandtl number). The appearance 
of this temperature ratio is associated with the variable property 
problem. But, the fact that only a temperature ratio appears is actually 
a considerable simplification; since for the general variable property 
fluid it is necessary to specify T w and T^ separately. 

Numerical solutions* of Eqs. (8a) for gas A were carried out for a 
wide range of values of T w /T_ for Pr = 0.7 and for selected values of 
T w /T„ for Pr = 1.0. The heat transfer results corresponding to these 
solutions are listed in Table II. The Nusselt and Grashof numbers, di- 
mensionless parameters, are given by Eqs. (12) and (13b), respectively. 


*A11 numerical integrations were carried out on an IBM Card Pro- 
grammed Calculator using a technique presented in detail in appendix B 
of reference 2. 


TABLE II 


Heat Transfer Results for Gas A 



, 1/4 U l/4 

Nu x,w/Gr x , w or Nu L,w/3 Gr L,w 

t 

O 

II 

& 

Pr = 1.0 

4 

0.371 " 


3 

0.368 

0.418 

5/2 

0.366: 


2 

0.363 


3/4 

0.348 


1/2 

0.339 


1/3 

0.330 

0.375 

1/4 

0.323:. 



Now, we proceed to generalize these results. Attention is first 
focused on Pr = 0.7. From an analytical solution for the constant 
property fluid , it is found that the heat transfer (for Pr = 0.7) is 

NUx/Gr*/ 4 = 0.353 = Nu L /| Gr^ 4 (19) 

Here, the Grashof number is given by Eq. (13a). We are immediately led 
to ask whether there is some way by. which the constant property result, 
Eq. (19), can be made to coincide with the variable property results 
appearing in Table II. It may be observed that the temperature for 
evaluating k, v, and B in Eq. (19) is at our disposal. With regard to 
P, it has already been shown that it is proper to replace p by l/T w 
when perfect gases are involved. It is a matter of making a few trials 

- This result may also be thought of as applying to gas A when “ 

V T - - 
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to determine the temperature for evaluating k and v of Eq. (19) which 
gives the best agreement between the constant property result and the 
variable property findings of Table II. This temperature, which may be 
termed a reference temperature T r , is found to be 

T r = T w - 0»38(T W - Tj, P = l/T — (20) 

As shown by the dashed line of Fig. 2, the error in the heat transfer 
predicted from the constant property result by using this reference tem- 
perature is at most 0.6 percent over the entire range l/4 < T w /T^ < 4. 
Figure 2 also shows the errors in the heat transfer predicted from the 
constant property result when k and y are evaluated at the wall tem- 
perature, at the ambient temperature, and at the film temperature T^.. 
p is taken as l/T„ in all cases. It is noteworthy that the use of T^, 
as a reference temperature provides heat transfer results which would be 
adequate for almost all engineering purposes. 

The results listed in Table II for Pr * 1 were obtained to check 
whether the reference temperatures found for Pr = 0.7 could be applied 
for the entire range of gas Prandtl numbers. When reference temperatures 
were computed for Pr = 1, they were found to be in excellent agreement 
with those for Pr = 0.7. So, it is felt that Eq. (20) is valid for the 
entire Prandtl number range of gases. 

Having established a reference temperature rule for gas A, we now 
proceed to see if it can be used for predicting variable property heat 
transfer results for other gases. 

Results for gas B . - As may be seen from Table II, gas B differs 
from gas A only by a lowering of the exponent in the conductivity and 
viscosity variations. The lower exponent becomes more appropriate in 
describing real gas behavior as the temperature increases. 
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Gas B will be treated in the following way: A special case will be 

selected and heat transfer results will be computed based on a solution 
of Eq. ( 8 a). Then, for the same special case, an alternate heat transfer 
result is computed for gas B using the constant property formula and the 
reference temperature rule of Eq. (20). The comparison of the two heat 
transfer results provides a test of the reference temperature procedure. 

The case selected for study is Pr = 0.7, T /T^ - 3. From a numeri- 
cal solution of Eq. ( 8 a), it is found that 

^ u x,w/^ r x{w = 0-373 =» Nu l >w / 3 Gr^v (21a) 

The Grashof number is given by Eq. (13b). Then, using the constant 
property formula, Eq. (19), and the reference temperature relation, 

Eq. (20), an alternate heat transfer result for gas B is 


Hu x,v/ Gr x{i - 0-STO = I»“L,w/t <*£(1 


( 21 b) 

It is seen that the prediction using the reference temperature procedure 
is very good. 

Results for gas C . - The very special mathematical simplifications 
which are associated with gas C have already been pointed outj see Eq. 
( 8 b) and associated discussion. Inspection of Eq. ( 8 b) shows that T w /T 
need not be specified. 

From a numerical solution of Eq. ( 8 b), the heat transfer result for 
gas C for Pr = 0.7 is found to be 


Nu x,w . Nu x,- . fVr _ 

Gr 1 / 4 Gr 1 / 4 Grp 4 

x,w x, — x,r 


= 0.353 


( 22 ) 
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where Gr x is given by Eq. (13b). A similar finding is true for the 
overall heat transfer. It is thus seen that, as a consequence of the 
special properties of gas C, the ratio Nu x /Gr x / 4 is independent of the 
second subscript indicating where the properties are evaluated. 

Looking now at the constant property formula, Eq. (19), it is seen 
that no matter how the properties are evaluated, its heat transfer result 
coincides with the variable property result of Eq. (22). (That is, pro- 
vided that (3 is evaluated as l/T^.) So, the reference temperature of 
Eq. (20) is certainly satisfactory for gas C. 

Results for gas D . - It has already been pointed out in the dis- 
cussion of Table I that a Sutherland type representation for k and p 
is generally closer to reality than are simple power laws. For gas D, 
the assumptions of constant Cp and Pr are maintained, so the same 
Sutherland relation is used for k and p. The constant A^^ appearing 
in Table I is assigned the value 362° F abs. This number is from a cor- 
relation by Glassman and Bonilla (4) of thermal conductivity data for air. 

With the aid of Table I, the property groupings pp and pk which 
appear in the differential Eqs. (8a) may be evaluated. When the dimen- 
sionless temperature 9 = (T - T^)/(T W - T ) is introduced, there results 



Careful inspection of this expression shows that T w and T must be 
separately specified before a solution of Eq. (8a) can be carried out. 
In other words, specific problems must be considered. 
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We proceed to test the reference temperature relation (20) in the 
same fashion described for gas B. Heat transfer results based on solu- 
tions of the differential equations are compared with alternate results 
determined using the reference temperature procedure. The following two 
special cases are considered^ 


Pr = 0.7 

r w = 1800° F abs V 


Pr = 

T = 600° F abs 


= Pr = 0.7 

case I T w = 600° F abs V case II 
T = 1800° F abs 


From numerical solution of Eq. (8a), the heat transfer results for these 


cases are 


(23a) 


1/4 — U 1/4 jO, 370, 

»“*,w/Gr x , u - Nu L)V /y G %w - { 


(23b) 


0.371, case I 
0.335, case II 

Alternately, heat transfer results can be computed for gas D for these 
cases by using the constant property formula, Eq. (19), and the reference 
temperature rule, Eq. (20). The results are 

case I 

0.335, case II 

Again, it is seen that the predictions using the reference temperature 
procedure are good indeed. 

Results for gas E . - As has already been pointed out, gas E is meant 
to be a close approximation to air. The k and p variations are rep- 
resented by Sutherland formulas (see table I) taken respectively from 
Glassman and Bonilla (4) and the BBS Tables (5). The constants Aj_ and 
Ag appearing in these representations are 362° F abs and 198,7° F abs 
(110.4° K) , respectively. The specific heat data, taken from the NBS 
Tables, is represented by two straight lines over the temperature range 
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studied. The variation of the Prandtl number need not he specified, 
since it is determined once Cp, k, and p are given. 

The. property groupings pn , pk, and c p/ c p which appear in the 
differential Eqs. (8a) may be evaluated with the aid of Table 1. Once 
this has been done, it is easy to see that it is necessary that T w 
and T„ be separately specified before proceeding witli a solution 
of Eqs. (8a). So again, as for gas D, specific problems must be 
considered. 

We now proceed to check the reference temperature procedure in the 
same manner as was already used for gases B, C, and D. The special 
cases selected here for this purpose are 


Pr w = 0.7 


T w = 1800 F abs 


T„ = 600° F abs J 


Pr w = 0.7 


Y Case I 


x w 


600° R 


Y Case II 


/ X / 4 — 
Nu x,w/ Gr x,w = Nu L 


U r x / 4 

>W / 3 Gr L,w = 


(24a) 


T. a 1800° R _ 

Numerical solutions of Eqs. (8a) provide the following heat transfer 
results 

* 

*0.358, case I 
0.346, case II 

We now proceed to compute alternate heat transfer results based on 
the constant property findings and the reference temperature relation, 
Eq. (20). In co nn ect! on with this computation, it is worthwhile noting 
a particular feature associated with the fact that the gas under con- 
sideration has a temperature- dependent Prandtl number. Over the tem- 
perature range studied for gas E, the Prandtl number variation was 4.5 
percent. First, it may be observed that the constant property analysis 
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yields 



function of Pr 


(25) 


where the Grashof number is given in Eq. (13a). This relation is plotted 
in Pig. 3 for the Prandtl number range appropriate to gases. Now, con- 
sider the application of Eq. (25) to a variable property fluid with 
temperature-dependent Prandtl number. Under these conditions, not only 
are k and v evaluated at the reference temperature, but also Pr. 

Using Eq. (25) and Fig. 3 in conjunction with the reference tem- 
perature relation (20), alternate heat transfer results for gas E for. - 
cases Itand II ate found to be 


l/4 / 4 l/4 I 0.361, 

^ u x,w/^ r x,w = ^ u L,w / 3 ® r L,w = I 


(24b) 


case I 

0,343 , case II 

Comparison of Eqs. (24a) and (24b) shows that even for this more com- 
plicated (and more realistic) situation, the reference temperature pro- 
cedure predicts very good heat transfer results. 


FLUID PROPERTIES OF LIQUID MERCURY 
As already noted in the INTRODUCTION, it was initially hoped to 
study the variable property problem in several liquid metals. However, 
it was observed that the property variations of the various liquid 
metals of technical interest were sufficiently dissimilar that separate 
investigations appeared to be necessary for each one. Moreover, the 
computing time required to obtain numerical solutions for liquid metals 
is an order of magnitude greater than that required for gases. So, it 
was found necessary to limit the current study to one liquid metal, 


namely, mercury. 
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When considering the property variations of liquid mercury, it might 
be well first to point out some qualitative trends and make comparisons 
with gases. For liquid mercury, the viscosity decreases with increasing 
temperature, while the thermal conductivity increases with temperature. 

On the other hand, for gases, both the viscosity and conductivity in- 
crease with temperature. The Prandtl number of mercury shows a strong 
(percentage) decrease with increasing temperature, while the Prandtl 
number of most gases varies by only a few percent over a large temper- 
ature range. Also, the absolute magnitude of mercury's Prandtl number- is 
1 or 2 percent of that of gases. The above remarks apply to normal 
engineering conditions, e.g., gas dissociation ’ is excluded. 

Property data for liquid mercury was taken from the Liquid Metals 
Handbook (6). To facilitate numerical integrations, polynomial repre- 
sentations of the following form were fitted to the data 


a n t n 

where t is in degrees Fahrenheit. The coefficients 
Table III. 



(26) 


an are given in 


* 


This is by no means the rule among liquid metals. 
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TABLE III 

Coefficients For The PolynomI cal Representations 
Of Properties Of Liquid Mercury 



k 


C P 

P 


Btu 

lb 

Btu 

lb 



ft - hr 

lb - °F 

cu ft 

*0 

4.47924 

4.34620 

0. 334620 (10' 1 ) 

851.514 

a l 

0. 830958 (10" 2 ) 

-0.991162(10“ 2 ) 

-0. 393353 (10' 5 ) 

-0.864880(l0 -1 ) 

a 2 

-0. 380163 (l0~ 5 ) 

0. 179060 (10 -4 ) 

0. 344649 (10' 8 ) 

0. 986194 (10“ 5 ) 

a 3 

0 

-0. 127524 (lO -7 ) 

0 

-0. 592566 (lO -8 ) 




HEAT TRANSFER RESULTS FOR LIQUID MERCURY 
First, based on numerical solutions of Eqs. (8), heat transfer re- 
sults are given for special cases. Then, these findings are utilized in 
determining reference temperatures. 

Results for specific cases . - For liquid mercury, the buoyancy force 
appearing in the first of Eqs. (8) does not simplify as it did for gases. 
So, it is necessary to deal with Eqs. (8) as they stand. Further, the 
nature of the property variations of liquid mercury requires that ty 
and t m be specified separately, i.e., that specific problems be con- 
sidered. The two cases chosen for study are described below. In both 
instances, t„ > t . 
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t w = 600° F 


t = 100° F 


Pr^ — 0 • 0085 


Pr = 0.022 


V case I 


^ = 450° F' 


t_ = 150 F 


Pr w = 0.010 


Pr_ = 0.019 


y case II 


Nu x,v/ Gr x{t = Nu L,w/ Gr i{t = 


(27) 


On the basis of the rather lengthy numerical solutions of the differ- 
ential Eqs. (8), the heat transfer results are found to bfee 

0.0501, case I 
^0.0556, case II 
where the Grashof number is given by Eq. (13). 

Reference temperature relation . - Now, it may be inquired as to 
whether there is some reference temperature rule which will cause the 
constant property heat transfer results to coincide with those of Eq. (27). 

It has already been noted that the constant property analysis yields 
heat transfer results which are given by Eq. (25). A plot of this re- 
lation is given on Fig. 4(a) for the Prandtl number range appropriate 
to liquid mercury. After making a few trials, it is found that by 
evaluating k, v, P, and Pr at 

T r = T w - 0.3 (T w - T.) (28) 

the constant property heat transfer results: of Fig. 4(a) coincide with 
those of Eq. (27). It is emphasized that this same reference temperature 
is found for both variable property cases studied here. 

It may next be inquired as to how well the film temperature, Tfj 
might serve as an alternate reference temperature. The heat transfer 
calculated for cases I and II by evaluating the constant property re- 
sults of Fig. 4(a) at the film temperature T f agree quite well with 
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Eq. (27), certainly within the range of most engineering requirements. 

The same statement applies when all the properties (including p) are 
evaluated at T r = T w - 0 ,38(T W - T^), which was the reference tempera- 
ture derived for gases. 

EXPERIMENTAL VERIFICATION 

A complete survey of available experimental heat transfer data for 
free convection on a vertical plate is presented in Ref. 1. It is noted 
there that the hulk of the experiments were carried out in air. 

For air, the conditions of most tests were such that 1< T w /T — <1.5. 
These experimenters took no note of variable property effects. This 
action is justified by the results of the present analysis; for example, 
see Fig. 2. In one experiment, that of Weise (7), a value of T w /T — of 
2.2 was achieved. Unfortunately, Weise 's apparatus, constructed for . 
horizontal plate tests, was not well suited for vertical plate 
studies. 

For liquid mercury, Saunders 1 (8) work represents the only experi- 
ment on a vertical plate. His tests were carried out at very small, 
temperature differences, and hence variable property effects did not 
enter. 

So, it appears that there are currently no data available to check 
the findings of the analysis presented here. 

RESULTS FOR BOUNDARY LAYER THICKNESS 
AND VELOCITY PARAMETERS 

While the heat transfer is by far the result of greatest practical 
importance, there are other quantities which may be of interest. A 
brief description is given below of several quantities for which results 
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are to be reported. Following the descriptive paragraphs, numerical 
findings derived from the constant property analysis are given, and then 
tables of reference temperature are supplied for extending these constant 
property results to variable property situations. 

(a) Boundary- layer thickness, 5^: The thickness of the boundary 

layer is by no means a precise concept, and its definition is somewhat 
arbitrary. The distance from the plate surface at which T - T has 
shrunk to a small fraction i of the overall temperature difference, * 

T w - T^, is used here to define a boundary layer thickness. This defi- 
nition, while satisfactory for Prandtl numbers below unity, should not 
be used for high Prandtl numbers. 

(b) Maximum velocity, u max : The vertical velocity .u (parallel 

to the plate surface) takes on zero values both at the surface and at 
y = •. and hence achieves some maximum between. 

(c) Location of the maximum velocity, y u : The distance from the 

plate surface at which the velocity maximum occurs is denoted by y u . 

(d) Friction coefficient, cp: The coefficient of friction provides 


a dimensionless presentation of the wall shear stress t w in the fol- 


lowing manner. 



( 29 ) 


The quotient, (v/x), plays the role of a characteristic velocity. 

(e) Flow rate, ¥: The upward flow generated by the free convec- 


tion forces is given by 



pu dy 


( 30 ) 


where b is the plate width. 
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Results from the constant property analysis . - Utilizing the tabu- 
lated constant property solutions of Ostrach (2) in conduction with those 
of Ref. 1, results for the quantities described in paragraphs (a) through 
(e) have been computed. A presentation of these findings for the Prandtl • 
number ranges appropriated to gases and to liquid mercury is made 
respectively on Figs. 3 and 4. Heat transfer results are also shown, The 
Grashof number is given by Eq. (l3a) . 

It was originally decided to define the boundary- layer thickness as 
the distance from the plate where (T - T„)/(T W - T„) = 0.02. For the 
Prandtl number range of gases, there was no difficulty evaluating such a 
definition from the available numerical solutions. But, for the low 
Prandtl number range, the numerical solutions were less precise at large 
distances from the wall, and it was necessary to use a thickness based on 
(T - T — )/'(T w - T — ) = 0.09. It is estimated that for the Prandtl number 
range of mercury, 8 q.02 * 1*25 &q.05* 

Variable property results; reference temperatures . - In a manner 
identical to that outlined for the heat transfer, the findings of the 
constant property analysis for the quantities described above may be ex- 
tended to the variable property situation. It (Jnly remains to present 
appropriate reference temperatures which are needed to evaluate the re- 
sults of Figs. 3 and 4. 

Such reference temperatures, computed using the variable property 
solutions which have been obtained here for gases and for liquid mercury, 
are listed in Tables if(^) and((l)). The heat transfer reference : temper- 
atures, already given by Eqs. ( 20) and (28), are included for completeness 
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TABLE IVU) 

Reference Temperature Relations For Gases 
(For use with Fig. 3) 


(T r - T w )/(T w - T w ) 


h 

0.38 

4 

i 

6 0.02 

0.67 

1 


^ax 

See below 

3 = 

-L 

2 T~ 

y u 

0.24 



c f 

0.10 



w 

0.85 

1 

f 


To find ^ax for gases, evaluate the Prandtl number of Fig. 3 at the 
mean of the extreme values of Pr in the particular problem under 
cons i deration . 

TABLE Iffy) 

Reference Temperature Relations 

For Liquid Mercury 
(For use with Figs. 4(a) and (b)) 

(T r - T W )/(T„ - T w ) 


h 

0.30 

8 0.05 

0.60 

^ax 

0 


0 

c f 

0.10 

w 

0.40 


CONCLUDING REMARKS 

From the findings reported here, it appears that free convection 
heat transfer under variable property conditions can be computed quickly 
and accurately by using the constant property results in conjunction with 
reference temperature relations. For gases and liquid mercury, the refer- 
ence temperature relations are given by Eqs. (20) and (28), respectively. 
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Further, it may be observed that the film temperature appears to serve 
as an adequate reference temperature (with 3 = l/T — for gases) for 
most engineering purposes. 

Reference temperature relations for use in computing boundary- layer 
thickness and velocity parameters for variable property conditions are 
given in Tables IV(a) and (b). 
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Figure 2. - Errors made in the prediction of the heat transfer for gas A by evaluating 
the constant property result at various temperatures. 
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Figure 3. - Results of the constant property analysis for the Prandtl number 
range of gases. 
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Figure 4(a) - Results of the constant property analysis for the Prandtl number 
range of liquid mercury. 
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